Ferromagnetic insulating state in tensile-strained LaCoO3_3 thin films

With local density approximation + Hubbard $U$ (LDA+$U$) calculations, we show that the ferromagnetic (FM) insulating state observed in tensile-strained LaCoO$_3$ epitaxial thin films is most likely a mixture of low-spin (LS) and high-spin (HS) Co, namely, a HS/LS mixture state. Compared with other FM states, including the intermediate-spin (IS) state (\textit{metallic} within LDA+$U$), which consists of IS Co only, and the insulating IS/LS mixture state, the HS/LS state is the most favorable one. The FM order in HS/LS state is stabilized via the superexchange interactions between adjacent LS and HS Co. We also show that Co spin state can be identified by measuring the electric field gradient (EFG) at Co nucleus via nuclear magnetic resonance (NMR) spectroscopy

Adaptive Genetic Algorithm for Crystal Structure Prediction

We present a genetic algorithm (GA) for structural search that combines the speed of structure exploration by classical potentials with the accuracy of density functional theory (DFT) calculations in an adaptive and iterative way. This strategy increases the efficiency of the DFT-based GA by several orders of magnitude. This gain allows considerable increase in size and complexity of systems that can be studied by first principles. The method's performance is illustrated by successful structure identifications of complex binary and ternary inter-metallic compounds with 36 and 54 atoms per cell, respectively. The discovery of a multi-TPa Mg-silicate phase with unit cell containing up to 56 atoms is also reported. Such phase is likely to be an essential component of terrestrial exoplanetary mantles.Comment: 14 pages, 4 figure

Metric tensor as the dynamical variable for variable cell-shape molecular dynamics

We propose a new variable cell-shape molecular dynamics algorithm where the dynamical variables associated with the cell are the six independent dot products between the vectors defining the cell instead of the nine cartesian components of those vectors. Our choice of the metric tensor as the dynamical variable automatically eliminates the cell orientation from the dynamics. Furthermore, choosing for the cell kinetic energy a simple scalar that is quadratic in the time derivatives of the metric tensor, makes the dynamics invariant with respect to the choice of the simulation cell edges. Choosing the densitary character of that scalar allows us to have a dynamics that obeys the virial theorem. We derive the equations of motion for the two conditions of constant external pressure and constant thermodynamic tension. We also show that using the metric as variable is convenient for structural optimization under those two conditions. We use simulations for Ar with Lennard-Jones parameters and for Si with forces and stresses calculated from first-principles of density functional theory to illustrate the applications of the method.Comment: 10 pages + 6 figures, Latex, to be published in Physical Review

Ensemble density-functional theory for ab-initio molecular dynamics of metals and finite-temperature insulators

A new method is presented for performing first-principles molecular-dynamics simulations of systems with variable occupancies. We adopt a matrix representation for the one-particle statistical operator Gamma, to introduce a ``projected'' free energy functional G that depends on the Kohn-Sham orbitals only and that is invariant under their unitary transformations. The Liouville equation [ Gamma , H ] = 0 is always satisfied, guaranteeing a very efficient and stable variational minimization algorithm that can be extended to non-conventional entropic formulations or fictitious thermal distributions.Comment: 5 pages, two-column style with 2 postscript figures embedded. Uses REVTEX and epsf macros. Also available at http://www.physics.rutgers.edu/~dhv/preprints/index.html#nm_meta

Evidence for a Peierls phase-transition in a three-dimensional multiple charge-density waves solid

The effect of dimensionality on materials properties has become strikingly evident with the recent discovery of graphene. Charge ordering phenomena can be induced in one dimension by periodic distortions of a material's crystal structure, termed Peierls ordering transition. Charge-density waves can also be induced in solids by strong Coulomb repulsion between carriers, and at the extreme limit, Wigner predicted that crystallization itself can be induced in an electrons gas in free space close to the absolute zero of temperature. Similar phenomena are observed also in higher dimensions, but the microscopic description of the corresponding phase transition is often controversial, and remains an open field of research for fundamental physics. Here, we photoinduce the melting of the charge ordering in a complex three-dimensional solid and monitor the consequent charge redistribution by probing the optical response over a broad spectral range with ultrashort laser pulses. Although the photoinduced electronic temperature far exceeds the critical value, the charge-density wave is preserved until the lattice is sufficiently distorted to induce the phase transition. Combining this result with it ab initio} electronic structure calculations, we identified the Peierls origin of multiple charge-density waves in a three-dimensional system for the first time.Comment: Accepted for publication in Proc. Natl. Acad. Sci. US

A New Scenario on the Metal-Insulator Transition in VO2

The metal-insulator transition in VO2 was investigated using the three-band Hubbard model, in which the degeneracy of the 3d orbitals, the on-site Coulomb and exchange interactions, and the effects of lattice distortion were considered. A new scenario on the phase transition is proposed, where the increase in energy level separation among the t_2g orbitals caused by the lattice distortion triggers an abrupt change in the electronic configuration in doubly occupied sites from an S=1 Hund's coupling state to a spin S=0 state with much larger energy, and this strongly suppresses the charge fluctuation. Although the material is expected to be a Mott-Hubbard insulator in the insulating phase, the metal-to-insulator transition is not caused by an increase in relative strength of the Coulomb interaction against the electron hopping as in the usual Mott transition, but by the level splitting among the t_2g orbitals against the on-site exchange interaction. The metal-insulator transition in Ti2O3 can also be explained by the same scenario. Such a large change in the 3d orbital occupation at the phase transition can be detected by linear dichroic V 2p x-ray absorption measurements.Comment: 5 pages, 5 figures, to be published in J. Phys. Soc. Jpn. Vol. 72 No. 1

Elasticity of Diamond at High Pressures and Temperatures

We combine density functional theory within the local density approximation, the quasiharmonic approximation, and vibrational density of states to calculate single crystal elastic constants, and bulk and shear moduli of diamond at simultaneous high pressures and temperatures in the ranges of 0-500 GPa and 0-4800 K. Comparison with experimental values at ambient pressure and high temperature shows an excellent agreement for the first time with our first-principles results validating our method. We show that the anisotropy factor of diamond increases to 40% at high pressures and becomes temperature independent.Comment: 10 pages, 3 figures, 1 tabl